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Abstract 

The  battle  trace  is  a  plot  representing  success  of  a 
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1.  Introduction 

The  "battle  trace"  is  a  high  level  indicator  of  success  in  a 
battle,  computed  as  a  function  of  time  into  the  battle.   Barr, 
Weir  and  Hoffman  [1991]  derive  the  battle  trace  based  on 
conditions  in  a  battle  assumed  to  follow  the  Lanchester  Square 
Law.   Suppose  combatants  X  and  Y  engage  in  a  battle  and  suppose 
the  battle  is  partitioned  into  portions  corresponding  to  time 
increments  Alf  A2,  A3,  .  .  .  ,  which  together  cover  the  duration  of 
the  battle.   In  the  ith  time  interval,  the  battle  trace  is 
defined  to  be  B=[ Ay/y] / [ Ax/x] ,  where  Ay  and  Ax  are  measures  of 
casualties  sustained  by  Y  and  X,  respectively,  within  the  given 
time  interval,  and  y  and  x  are  measures  of  the  strengths  of  Y  and 
X,  respectively,  during  that  period. 

Barr,  Weir  and  Hoffman  [1991]  also  consider  two  symmetric 
variants  of  B,  also  refered  to  as  "battle  trace  measures:" 

the  "symmetric  battle  trace," 

f  B,  if  B>1; 
Bs=  { 

{   2-1/B,  if  B<1, 

and  the  "log  trace," 

BL=ln(Ay+a)  +  ln(x+a)  -  ln(y+a)  -  ln(Ax+a), 

where  a  is  a  suitable  positive  constant.   (Battle  traces  with 

arguments  increased  by  such  a  are  called  "incremented"  in  what 

follows.)   They  discuss  how  these  measures  might  be  used  in 

applications  related  to  field  testing  and  combat  simulation 

models. 
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In  the  present  paper  we  discuss  possible  uses  of  battle 
traces  in  combat  training  environments.   Specifically,  we 
envision  that  the  Janus  combat  simulation  model  might  be  used  to 
enhance  training  conducted  at  the  National  Training  Center  (NTC) 
at  Fort  Irwin,  CA.   Once  a  unit  has  completed  a  set  of  battles  at 
the  NTC,  the  relevant  parameters  of  the  battles  (such  as 
scenarios,  routes  followed  by  forces  on  each  side,  uses  of 
scouting  platoons,  engineering  activities,  etc.)  can  be  set  (with 
varying  degrees  of  fidelity)  into  the  Janus  combat  simulation, 
and  the  battles  can  be  replayed  within  Janus.   The  objective  of 
such  replays  is  to  provide  the  unit  commander  information  about 
how  the  battle  might  have  progressed  had  various  changes  been 
made  by  the  commander  in  his  prosecution  of  the  battle.   For 
example,  a  commander  might  wish  to  examine  how  his  forces  would 
have  done  had  he  used  his  scouting  assets  more  effectively,  or  if 
he  had  exploited  terrain  masking  more  effectively  in  maneuvering 
his  forces  toward  an  attack. 

However,  there  is  considerable  uncontrolled  variation  in 
battle  outcomes,  both  for  battles  at  the  NTC  and  those  simulated 
within  Janus.   Thus,  even  if  one  replicates  a  battle  by  holding 
the  controllable  factors  at  the  same  levels  while  running  a 
second  battle,  generally  there  will  be  substantial  differences 
between  the  two  battles  nonetheless.   Any  attempt  to  exactly 
reproduce  the  results  of  a  given  battle  is  doomed  to  failure,  so 
by  extension  there  is  not  a  one-to-one  relationship  between  the 
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set  of  values  of  parameters  in  a  battle  simulation  and  consequent 
ways  in  which  the  simulated  battles  unfold  over  time.   It  may  be 
difficult,  therefore,  for  a  battle  commander  to  interpret  effects 
on  battle  results  that  might  be  due  to  changes  in  battle 
parameters  or  conditions.1 

With  the  substantial  variations  in  battle  outcomes,  how  can 
a  commander  get  a  clear  picture  of  the  likely  effects  of  changing 
selected  battle  parameters?   Estimates  based  on  traditional 
measures  of  effectiveness  (MOE's)  generally  provide  definitive 
information  about  battle  outcomes  only  when  it  is  feasible  to 
carry  out  numerous  replications  with  each  combination  of  factors 
under  study.   We  doubt  this  is  feasible  in  the  context  of  using 
Janus  to  augment  NTC  training  as  mentioned  above.   Thus  it  is 
desirable  to  employ  alternative  measures  that  give  general, 
"high-level"  views  of  the  progress  of  each  battle  while 
atttempting  to  avoid  the  volatility  associated  with  specific  end- 
of-battle  MOE's.   We  believe  battle  traces  may  be  useful  in  this 
regard. 

Several  questions  about  computation  of  battle  trace  measures 
were  investigated: 


1This  is  precisely  why  f orce-on-f orce  field  tests  generally 
are  designed  to  employ  many  replications  within  each  cell  of  the 
test  matrix.  Such  "variance  reduction  by  replication"  is  not 
feasible  in  most  training  applications  at  the  present  time  because 
of  resource  constraints.  Eventually  it  may  be  possible  to  run  a 
significant  number  of  Janus  replications  simultaneously  using 
parallel  computing  techniques,  so  meaningful  sample  sizes  can  be 
attained  in  "real  time." 
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*  How  should  the  time  increments  At,  be  chosen? 

*  Which  version  of  the  battle  trace  is  best  for  training 
applications? 

*  Should  the  battle  trace  plot  be  smoothed?   If  so,  what  is 
the  best  method  of  smoothing? 

Numerous  crude  simulation  experiments  were  conducted  to 
provide  rough  initial  answers  to  these  questions.   These 
experiments  provided  information  relevant  to  making  a  choice  of 
which  battle  trace  measure  should  be  used  in  the  initial 
demonstrations  of  the  concept.   The  resulting  choice  is  described 
below.   Whether  this  particular  measure  does  indeed  help  a 
commander  understand  the  general  effects  of  making  certain 
changes  in  how  he  fights  a  given  battle  is,  at  this  point,  an 
open  question.   This  issue  can  be  examined  by  future 
experimentation.   An  indicator  of  the  utility  of  the  measure  will 
be  the  degree  to  which  it  is  readily  accepted  and  used. 

In  the  sections  that  follow,  we  describe  properties  of  the 
selected  form  of  the  battle  trace  and  document  how  it  was 
implemented  in  the  Janus (A)  simulation  at  TRAC-Monterey .   We 
expect  demonstrations  with  this  form  of  the  battle  trace  may  lead 
to  performance  of  experiments  to  evaluate  its  utility  in  the 
training  context. 
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2.  Initial  Simulation  Results 

As  mentioned  above,  early  in  our  investiagtion  we  generated 
various  battle  traces  with  data  obtained  from  crude  simulations. 
The  simulated  battles  followed  the  Lanchester  Square  Law  within 
each  time  interval,  Atj,*  i=l,2,...,n,  in  a  partition  of  the 
battle  period.   This  was  done  to  generate  data  required  in 
computation  of  the  battle  trace  in  its  various  forms.   The 
general  approach  was  to  examine  the  appearance  and  stability  of 
various  forms  of  battle  trace  displays  computed  with  the 
generated  data.   Several  related  mathematical  developments  are 
presented  in  Appendices  1-3.   In  Section  6  we  report  results  of 
preliminary  experiments  with  the  selected  battle  trace  when  it  is 
run  from  within  Janus. 

We  drew  the  following  general  conclusions  from  these  initial 
experiments : 

(1)  The  plots  of  the  various  battle  trace  functions  have 
considerable  variation  (i.e.,  are  not  smooth).   This  is 
aggravated  when  the  time  sub-intervals  At  are  too 
small.   We  found  that,  for  the  simulated  battles  we 
considered,  plots  with  20  to  50  data  points  for  the 
entire  battle  (lasting  about  40  minutes)  looked  best. 

(2)  Potential  division  by  zero  in  B  and  Bs,  and  zero 
arguments  within  log  functions  in  BL,  cause  undesirable 
and  unnecessary  problems  in  the  computations.   It 
appears  the  "incremented"  versions  of  the  battle  trace 
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functions  successfully  avoid  these  problems  without 
significant  loss  of  fidelity. 
(3)   Plots  using  dynamic  definitions  of  At  appear  to  have 

advantages  over  those  with  periods  of  fixed  length.  It 
appears  definitions  involving  Ax  or  Ay  being 
sufficiently  large  (or  attaining  sufficiently  large 
values  of  other  quantities  related  to  expected 
casualties,  such  as  the  number  of  shots  fired  or 
received  by  one  or  both  forces)  work  well.   The  values 
of  thresholds  defining  "sufficiently  large"  should  be 
chosen  so  as  to  provide  an  incremented  battle  trace 
with  a  fairly  smooth  plot.   Choice  of  such  a  threshold 
generally  involves  the  expected  length  and  outcome  of 
the  battle  and  the  numbers  of  combatants  involved.  It 
appears  that  not  more  than  about  50  data  points  should 
be  generated  over  the  duration  of  the  battle;  fewer 
data  points  would  be  appropriate  if  there  are 
relatively  few  combatants  involved. 
As  a  result  of  our  initial  experiments  we  concluded  that,  on 
balance,  an  incremented  version  of  the  log  trace,  BL  is  a  good 
choice.   Points  on  the  log  trace  should  be  computed  with  At 
defined  dynamically.   For  runs  conducted  within  Janus (A)  at  TRAC- 
Monterey  we  implemented  a  rule  involving  the  number  of  shots 
fired  by  or  received  by  each  task  force  under  consideration.   The 
purpose  of  computing  battle  trace  in  Janus (A)  was  to  conduct 
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demonstrations  of  the  concept  and  to  provide  further  experiments 
and  evaluations  of  the  utility  of  battle  traces.   Details  of  this 
implementation  are  given  in  the  Section  6. 

3.  Discussion  of  the  Recommended  Battle  Trace  Form 

Even  though  we  are  recommending  the  log  trace  BL,  we  discuss 
computation  of  B  here,  since  its  form  admits  an  interesing  and 
relatively  simple  interpretation.   Since  BL  is  simply  related  to 
B,  this  discussion  is  relevant  to  BL  as  well.   We  will  use  the 
term  "battle  trace"  in  a  generic  sense  to  refer  to  any  form  of 
the  measures  under  discussion. 

The  numbers  of  casualties  occuring  on  each  side  during  an 
interval  of  time  can  be  considered  to  be  outcomes  on  random 
variables.   The  means  of  these  random  variables,  that  is,  the 
expected  numbers  of  casualties  in  the  time  interval,  can  be 
estimated  by  summing  the  pk '  s  associated  with  shots  upon  each 
force  during  the  time  interval.   Such  estimators,  called  "sums  of 
pk's"  and  denoted  "Zpk,"  are  discussed  in  Appendix  2.   There,  it 
is  shown  that  use  of  appropriate  Zpk  •  s  to  represent  the  expected 
casualties  Ax  and  Ay,  will  reduce  the  variance  of  computed  values 
of  points  on  battle  traces,  and  hence  will  make  more  precise  (and 
possibly  appear  to  "smooth")  their  plots.   We  used  this  technique 
in  forming  the  log  trace  discussed  below. 

In  a  similar  way,  the  numbers  of  relevant  combatants 
involved  in  the  battle  at  a  given  time  (i.e.,  values  of  y  and  x 
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in  the  formula  for  BL)  can  be  considered  to  be  outcomes  on  random 
quantities.   It  is  desirable  to  use  variance  reducing  estimators 
for  the  means  of  the  distributions  of  x  and  y  within  each 
corresponding  battle  time  interval.   By  analogy  with  the 
discussion  in  Appendix  2,  it  appears  that  an  estimator  similar  to 
the  sum  of  pk ' s  estimator  would  be  appropriate  for  this  purpose. 
Consider  estimating  a  quantity  to  represent  x  in  the  battle  trace 
computation.   Within  the  Lanchester  Square  Law,  x  denotes  the 
size  of  the  X  force  at  the  time  in  question;  note,  however,  the 
effect  of  these  X  forces  is  to  decrease  y  by  approximately  b-x. 
The  effect  of  x  in  the  Lanchester  model  is  thus  related  to  the 
number  of  units  in  the  X  force  that  may  engage  Y  units.   One 
possibility  for  the  battle  trace  computation  is  thus  to  count,  in 
the  ith  time  interval,  the  number  of  X  units  that  have  line  of 
sight  (or,  alternatively,  have  detected  or  have  acquired)  at 
least  one  unit  of  the  Y  force.   This  would  correspond  roughly  to 
the  body  count  estimator  (as  discussed  in  Appendix  2)  of  the  mean 
of  the  x  distribution. 

As  an  alternative,  we  propose  summing  an  appropriate  measure 
of  the  potential  of  each  X  unit  to  inflict  casualties  on  the  Y 
force.   Due  to  the  way  in  which  the  Janus  simulation  processes 
information  and  updates  target  lists,  it  is  feasible  to 
determine,  for  each  unit  of  the  X  force,  the  maximum  probability 
of  kill,  maxp,  within  the  ith  time  interval,  against  any  target 
which  is  detectable  by  the  unit  (i.e.,  is  intervisible,  is  within 
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the  unit's  sector  of  regard,  and  is  within  its  maximum  range 
limits)  .   Let  maxpj  denote  the  value  of  maxp  for  the  jth  X-unit  in 
the  ith  time  interval  and  suppose  there  are  N  X-units  alive  at 
the  beginning  of  this  time  interval.   We  propose  using  the  value 
EJ=1Nmaxpj  to  estimate  the  mean  of  k-  x,  and  to  use  this  in 
computation  of  the  battle  trace.   Note  this  value  is  a  measure  of 
the  potential  of  the  X  force  to  inflict  casualties  during  the  ith 
time  interval.   A  similar  estimator,  summing  maxp  over  Y-units, 
measures  the  potential  of  the  Y  force  to  inflict  casualties  on  X. 
This  is  used  to  estimate  the  mean  of  k-y.   In  the  battle  trace 
computation,  the  unknown  proportionality  constant  k,  assumed  to 
be  the  same  in  numerator  and  denominator  of  B,  cancels. 

If  the  sum  of  pk '  s  estimators  are  used  for  Ax  and  Ay  and  the 
sum  of  maxp's  are  used  for  x  and  y  as  described  above,  the  battle 
trace  can  be  interpreted  in  the  following  way.   Let  pxk  denote 
the  probability  of  kill  by  an  X-unit  (against  a  Y-unit)  in  the 
kth  engagement,  with  a  similar  definition  for  pyk;  let  maxpXj 
denote  the  maximum  probability  of  kill  by  the  jth  X-unit  against 
any  of  its  Y-targets  during  the  ith  time  interval,  with  a  similar 
definition  for  maxpyj .   Then  for  the  ith  time  interval,  the  battle 
trace 

B=[Ay/y]/[Ax/x] 
is  computed  as 
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£     p*  £     py 

p    „        shots  on  Y /        shots  on  X 


£     maxpy  £     maxpx 

Y-  units  X- units 


£       px  £     maxpx 


shots  on  Y  X-uni  ts ( f,  \ 


£       PY  £     maxpy 

shots  on  X  Y-units 

The  first  of  the  ratios  in  equation  (6)  is  the  expected 
casualty  exchange  ratio,  and  is  a  measure  of  the  relative  firing 
effectiveness  of  the  X  and  Y  forces  during  the  ith  time  interval. 
The  numerator  (denominator,  respectively)  is  a  sum  of  single  shot 
probabilities  of  kill  for  all  shots  taken  by  the  X  force  at  the  Y 
force  (all  shots  by  the  Y  force  on  the  X  force)  in  the  given  time 
interval.   These  terms  account  for  the  influence  of  weapon  firing 
during  an  interval  of  battle  which  includes  both  direct  and 
indirect  (artillery)  round  impacts.   Because  there  are  typically 
many  more  shots  fired  than  kills,  such  values  contain  more 
information  about  the  related  battle  processes  than  does  a  simple 
count  of  casualties. 

The  second  of  the  ratios  in  equation  (6)  is  a  force 
capabilities  ratio;  it  measures  the  relative  attrition  potential 
of  the  X  and  Y  forces.   The  quantity  in  numerator  of  this  ratio 
is  the  sum  of  the  maxumum  values  of  single  shot  probability  of 
kill  for  each  X  weapon  system  which  can  detect  at  least  one 
system  of  the  opposing  (Y)  force  some  time  during  the  given  time 
interval.   Conditions  on  detectability  include  existence  of  line 
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of  sight  between  the  combatants,  as  well  as  conformance  with 
bounds  on  the  sector  of  regard  (the  X-unit's  search  sector)  and 
range. 

Similar  comments  hold  for  the  denominator  of  the  second 
ratio  in  (6).   Thus  these  sums  contain  information  on  the 
potential  of  a  force  to  influence  the  outcome  of  battle 
irrespective  of  weapon  firing  events  and  independent  of  rate  of 
fire.   These  values  may  be  seen  to  reflect  the  mass  of  force  on 
each  side  of  the  conflict  generated  within  each  time  interval  of 
the  battle.   A  weapon  system  which  has  run  out  of  ammunition  or 
has  no  ammunition  suitable  for  engaging  detectable  targets  can 
make  little  direct  contribution  to  the  battle.   This  fact  is 
reflected  in  the  second  ratio  in  equation  (6)  since  then  the 
corresponding  maxp  is  zero,  hence  such  a  system  is  not  counted  in 
the  expressions  corresponding  to  x  or  y  in  the  formula  for  B. 
Selecting  the  maximum  single  shot  probability  of  kill  value  as 
the  potential  contribution  for  a  given  weapon  system  is 
consistent  with  an  assumption  that  each  weapon's  crew  will  work 
to  maximize  its  contribution  to  the  battle. 

In  terms  of  the  two  ratios  in  equation  (6)  above,  the  X 
force  is  winning  the  battle  at  a  given  point  in  the  battle  if  it 
is  engaging  targets  more  effectively  than  is  Y,  and  it  is 
maneuvering,  positioning,  using  cover  and  concealment  and  in 
general  is  preparing  for  battle  more  effectively  than  is  its 
adversary.   These  conditions  will  generate  relatively  high  values 
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of  the  casualty  exchange  ratio  (first  part  of  equation  6)  and  the 
attrition  potential  ratio  (second  part  of  equation  6) ,  so  the 
battle  trace  B  will  be  significantly  larger  than  1.0.   If  the 
incremented  log  trace  is  being  used  (as  we  recommend) ,  BL  will 
then  be  significantly  larger  than  zero. 

On  the  other  hand,  X  may  be  losing  at  a  given  point  in  the 
battle  even  though  it  enjoys  a  superior  casualty  exchange  at  that 
point  (first  ratio  in  (6)  greater  than  1.0),  if  Y  is  out- 
maneuvering  and  out-planning  X  so  as  to  have  sufficiently  high 
advantage  in  attrition  potential  (second  ratio  in  (6)  much 
smaller  than  1.0) .   It  is  then  possible  that  the  product  B  of  the 
two  ratios  is  less  than  1.0,  so  the  log  trace  would  be  less  than 
zero. 

4.  Implementation  Within  Janus (A)  at  TRAC-Monterey 

In  Janus  simulations  of  battle,  forces  on  each  side  may  be 
partitioned  into  task  forces  by  the  modeler.   The  compositions  of 
these  task  forces  can  be  varied  by  a  force  commander  during  the 
play  of  a  battle.   We  next  consider  computation  of  a  battle  trace 
for  each  of  possibly  several  task  forces  on  each  side,  where  the 
task  force  compositions  at  each  point  of  the  battle  are 
consistent  with  any  changes  made  by  the  force  commanders  during 
Janus  execution.   To  simplify  our  discussion,  first  we  consider 
computation  of  the  battle  trace  for  the  jth  Blue  Task  Force,  Bj , 
which  we  consider  to  be  on  side  X.   Let  us  begin  by  describing 
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how  each  of  the  sums  in  equation  (6)  are  computed  for  this  task 
force. 

The  sum  of  pk '  s  corresponding  to  the  first  term  in  the 
battle  trace,  Ay,  is 


Pk(B-R) 


all  red  who  can  see  B 

but  are  shot  at  by 

any  Blue  unit 


where  "Pk(B,R)"  denotes  a  probability  of  kill  by  a  Blue  firer  upon 
a  Red  target.   Similarly,  the  sum  of  pk '  s  corresponding  to  the 
third  term  of  B,  Ax,  is 


2_,    Pk(R-Bj) 
all  units 
€  B, 

The  sum  of  maxp/s  corresponding  to  the  second  term  of  B,  yf 


is 


maXPk{R-3) 


all  Red  uni ts 
for  which  at  least  one 


and  the  sum  of  maxpk's  corresponding  to  the  fourth  term  of  B,  x, 
is 

£        maXPk(B-R) 
units  e  B:  who 
could  detect  any 
red  uni t 

A  separate  battle  trace  should  be  computed  for  each  task 
force  of  interest  on  the  Blue  side  and  each  task  force  of 
interest  on  the  Red  side.   Additionally,  combinations  of  task 
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forces  could  be  considered  (all  Blue  units,  for  example)  or 
pseudo  task  forces  might  be  defined  (all  Blue  tanks,  for 
example) .   In  the  following  section  we  describe  how  the 
quantities  in  equation  (6)  were  captured  within  the  Janus (A) 
simulation  for  each  task  force  on  each  side. 

5.  Details  of  Implementation 

This  section  contains  a  description  of  how  the  Janus (A)  code 
was  modified  at  TRAC-Monterey  so  experiments  and  demonstrations 
with  the  battle  trace  could  be  carried  out.   Specific  listings  of 
Fortran  code  comprising  major  subroutines  of  Janus (A)  that  were 
so  modified  are  shown  in  Appendix  5.   Together  with  the  comments 
in  this  section,  the  listings  in  Appendix  5  should  serve  as 
documentation  of  our  implementation  of  battle  trace  in  Janus. 

With  the  revised  Janus  code,  several  arrays  are  maintained 
throughout  a  simulation  of  a  battle,  to  provide  the  data  required 
to  support  computation  of  battle  traces  at  the  end  of  appropriate 
time  intervals.   These  are  defined  as  follows,  where  "numsides" 
is  the  number  of  sides  (l=Blue  and  2=Red) ,  "numtasks"  is  the 
number  of  task  forces  on  a  given  side  and  "numunits"  is  the 
number  of  units  on  a  given  side:2 

SSHOTPKS(nurfeides)x(numtasks)  -  the  ijth  element  is  the  sum  of 

pk ' s  of  all  shots  by  side  i  against  all  elements  of  side  3-i 


The  dimensions  of  each  array  are  shown  as  subscrpts  to  the 
array  name. 
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who  could  detect  at  least  one  unit  in  task  force  j  of  side 
i.   The  elements  of  this  array  are  updated  just  before  each 
call  to  subroutine  SHOT  within  subroutine  RELOAD  of  Janus. 
The  update  involves  adding  the  pk  value  to  be  used  in  SHOT 
to  certain  components  of  the  row  of  SSHOTPKS  corresponding 
to  the  side  of  the  firing  unit.   The  components  that  are 
incremented  are  those  corresponding  to  side  i  task  forces 
detectable  by  the  target  fired  upon. 

SIPCTPKS(numsides)x(nuTltasks)  -  the  ijth  element  is  the  sum  of 
pk ' s  of  shots  fired  at  the  task  force  j  of  side  i  (by  any 
f irer  on  side  3-i) .   The  elements  of  this  array  are  updated 
just  prior  to  each  shot  by  adding  the  pk  for  the  shot  called 
in  subroutine  SHOT,  as  described  above,  to  the  appropriate 
element. 

BTRACEVAL(nunisides)x(numunJts)x(numtasks)  -  the  ijkth  element  is  the 
maximum  of  all  pk  values  for  the  jth  unit  on  the  ith  side 
against  any  detectable  units  of  the  kth  task  force  of  the 
(3-i)th  (opposing)  side.   This  array  is  updated  each  time  a 
target  list  is  determined  within  Janus  (in  subroutine  DETECT 
of  Janus) . 

BTINFL(numsides)x(numunJts)x(nuntasks)  -  the  ijkth  element  is  a  flag 
indicating  whether  the  jth  unit  on  side  i  could  detect  the 
kth  task  force  on  the  (3-i)th  side.   This  array  is  updated 
each  time  a  target  list  is  determined,  as  described  above. 

KNUMSHOTS(ni™,des)x(nuTltasks)  "   th«  ijth  element  is  the  number 
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of  shots  fired  by  or  on  the  jth  task  force  of  side  i. 

BTMAX(numsides)x(nlmJnits)  "   the  ijth  element  is  the  maximum 
over  columns  (opposing  task  forces)  of  BTRACEVAL  (i.e.,  the 
max-reduction  of  BTRACEVAL  across  columns) .   This  reduction 
is  calculated  just  prior  to  computing  a  point  on  the  battle 
trace  for  the  jth  task  force  of  side  i. 

PKUNITMAX(numsJdes)x(numtasks)  -  the  ijth  element  is  the  sum  of 
elements  in  the  jth  column  of  BTRACEVAL,  giving  the  sum  over 
units  on  side  3-i  of  max  pks  against  task  force  j  of  side  i. 
Elements  of  this  array  are  updated  when  BTMAX  is  updated. 

PKTASKMAX(nuT,s.des)x(nunitasks)  -  the  ijth  element  is  the  sum  of 
BTMAX(i,k)  over  k  (units)  in  taskforce  j. 

The  battle  traces  are  computed  using  values  in  the  arrays 
described  above.   At  each  cycle  of  Janus  (every  20  seconds  of 
battle  time) ,  computation  of  a  point  on  the  battle  trace  for  task 
force  JTASK  on  side  JSIDE  is  undertaken  if  the  condition 
KNUMSHOTS(JSIDE, JTASK)  >  t,  where  the  threshold  t  is  set  as  an 
input  parameter  in  the  Janus  set-up.   Use  of  this  rule  provides 
dynamic  allocation  of  the  lengths  of  the  time  intervals  over 
which  the  battle  traces  are  computed;  the  times  will  occur  as 
discrete  multiples  of  the  20  second  cycle  time  inherent  in  the 
Janus  update  schedule.   Points  will  be  computed  more  frequently 
for  a  task  force  when  it  is  heavily  engaged  with  the  enemy 
(either  firing  or  receiving  fires) . 

When  threshold  conditions  are  met  for  one  or  more  task 
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forces,  the  corresponding  points  on  the  battle  trace  plots  are 
computed  and  re-initialization  of  appropriate  elements  of  the 
arrays  described  above  takes  place. 

When  a  threshold  condition  is  met,  BTMAX  and  PKUNITMAX  are 
updated  and  the  components  needed  for  computation  of  the 
corresponding  battle  trace  point  are  extracted  from  the  arrays 
described  above.   Following  the  battle  trace  computation 
appropriate  elements  of  the  arrays  are  set  to  zero  in  preparation 
for  accumulating  values  in  the  next  time  interval  for  the 
succeeding  point  on  that  battle  trace.   In  terms  of  the 
quantities  in  the  arrays  described  above,  the  battle  trace  value 
shown  in  equation  (6),  for  task  force  j  of  side  i,  becomes: 

B   .  SSHOTPKS(i,j)  PKTASKMAX ( i , j) 

SIPCTPKS(i,j)     X  PKUNITMAX (i , j) 

Immediately  after  computing  the  point  on  the  battle  trace,  the 
following  array  elements  are  set  to  zero: 

SSHOTPKS(i, j)  ; 

SIPCTPKS(i, j)  ; 

BTMAX ( i , j )  ; 

BTRACEVAL(i, j,k)  for  all  task  forces  k  on  side  3-i. 

The  following  pseudo-code  illustrates  how  this  is 
implemented;  specific  listings  of  the  fortran  code  are  given  in 
Appendix  5. 
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Pseudo-code  for  computing  log  trace 

**comment:  check  each  battle  trace  threshold  condition** 

for  i=l  to  numsides 

for  j=l  to  numtasks 

if  KNUMSHOTS ( i , j )  <  t  then  go  to  GOON 

**comment:  if  condition  met,  compute  BL  ** 

R=0 

for  k=l  to  numunits      **comment : sum  column  of 

BTRACEVAL  &  initialize** 
R=R+BTRACEVAL ( 3 -i , k , j ) 
BTRACEVAL ( 3 -i , k , j ) =0 
next  k 
B=0 

for  L=l  to  numunits    **comment : sum  maxpk  for 

units  in  taskforce  & 
initialize** 
if  L  in  taskforce (i, j)  then  B=B+MAXPK ( i , L) 
if  L  in  taskforce (i, j)  then  MAXPK ( i , L) =0 
next  L 
BL(i, j)=ln(SSHOTPKS(i, j)  +  1)  -  In (SIPCTPKS ( i , j )  +  1) 

-ln(R  +  1)  +  ln(B  +  1) 
SSHOTKPKS(i, j)  =  0 
SIPCTPKS (i,j)  =  0 
GOON  next  j 
next  i 


6.  Results  of  Preliminary  Experiments  with  Janus (A) 

To  gain  a  tentative  impression  of  battle  traces  obtained 
with  the  modified  Janus (A)  simulation,  we  conducted  limited 
experiments  at  TRAC-Monterey .   These  experiments  all  relate  to 
battles  simulated  with  a  "Fulda  Gap"  scenario  developed  at  TRAC- 
Monterey  for  use  with  a  Janus  tutorial.   We  obtained  intuitive 
confirmation  of  the  indications  gained  with  the  earlier  crude 
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simulation  study  described  in  Section  2:  the  incremented  log 
trace  with  dynamic  definition  of  At  appears  to  be  a  good 
candidate  for  more  extensive  testing  in  a  training  context. 

We  observed  that  there  may  be  considerable  variation  in 
specific  points  on  battle  traces  for  iterations  of  the  same 
battle.   However,  the  general  features  of  battle  traces  for  such 
iterations  appear  to  be  quite  similar,  suggesting  the  battle 
traces  are  capable  of  capturing  major  events  and  trends  of  a 
battle.   Plots  of  battle  traces  for  iterations  of  the  Fulda  Gap 
scenario  are  shown  in  Appendix  4.   We  believe  the  variations  seen 
in  battle  traces  for  iterations  of  a  battle  simply  show  how 
different  the  specific  details  of  battles  can  be,  even  though  the 
controllable  parameters  are  held  constant  (see  Figure  3  of 
Appendix  4) . 

It  was  observed  that  the  log  trace  tended  to  have  a  shape 
similar  to  a  plot  of  casualty  exchange  ratios  (CER)  calculated 
within  each  time  interval  (Figure  5) .   This  is  to  be  expected 
since  the  CER  is  one  of  two  ratios  involved  in  computation  of  the 
battle  trace  (see  equation  (6)).   But  there  is  theoretical 
support  for  this  relationship,  based  on  the  Lanchester  model. 
This  can  be  seen  as  follows:  assuming 

dB/dt  =  -BR   and  dR/dt  =  -aB 
it  follows  that 
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Thus 


-^  /  -a 
_B  dt  =_£XA* 

R  dB    I  _o       a      AB 

dt 


Battle  trace  -  4^*-   s  -^  *  <  4^  ) 2 
AS     R         a         AB 

so 

B .   =  In  (-2.)    +  2ln(-^)    =  C0  +  c1ln(CB^)      . 
^  a  AB  ° 

If  the  Lanchester  model  fits  reasonable  well,  a  plot  of  log  trace 
versus  ln(CER)  should  appear  as  a  line  with  slope  c1  =  2 .   If 
this  appears  tenable,  one  can  estimate  the  ratio  (a/B)  by  eco. 
We  show  a  plot  of  BL  versus  ln(CER)  in  Figure  6,  together  with  a 
line  with  slope  2.   It  appears  the  line  fits  the  data  from  the 
Fulda  Gap  scenario  rather  well. 

Residuals  representing  differences  between  plotted  points  of 
the  form  (BL,ln(CER))  from  the  line  with  slope  2  show  variations 
in  (B/a)  at  the  end  of  each  time  interval  of  battle.   That  is, 
variations  of  the  plotted  points  from  the  line  indicate 
variations  of  the  battle  away  from  a  Lanchester  square-law 
battle.   Indeed,  one  can  get  rough  estimates  of  ln(B/a),  within 
time  increments  making  up  the  battle  period,  using  the  values 
(residual  +  c0)  ,  where  c0  is  obtained  from  the  fit  of  the  line 
with  slope  2.   A  plot  of  these  values  versus  time  might  serve  as 
an  indicator  of  how  well  Red  is  doing  relative  to  Blue;  we  have 
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plotted  such  curves  in  Figure  7  for  seven  iterations  of  the  Fulda 
Gap  scenario. 

In  a  masters  thesis  written  at  the  Air  Force  Institute  of 
Technology,  H.K.  Choi  suggested  a  variant  form  of  the  battle 
trace  which  he  thought  might  have  greater  stability  than  has  BL. 
Mr.  Choi  suggested  that  since  Blue  wins  in  a  Lanchester  square 
law  battle  whenever  C  =  BR2  -  aB2  <  0,  it  might  be  useful  to  plot 
values  of 

-dB  -dR 


C  -   -J±-xR*   _  _J*t_xB*  „  ***  -  AR-B 
R  B  At 

within  each  time  interval  versus  time  into  the  battle.   We  think 
this  suggestion  has  potential  merit,  and  we  have  plotted  this 
form  of  battle  trace,  which  we  denote  "Bc,"  for  seven  iterations 
of  the  Fulda  Gap  scenario,  in  Figure  9. 
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Appendix  1.  A  Time  Transformation  in  the  Lanchester  Square  Law 

In  the  initial  crude  simulations  used  to  generate  data  for 
the  battle  trace  experiments  described  in  Section  2,  we  varied 
only  one  attrition  parameter,  a(t),  in  the  Lanchester  Square  Law, 
The  following  argument  shows  this  can  be  done  without  loss  of 
generality. 

Consider  a  transformation  T=f(t)  of  the  time  parameter  in 
the  Lanchester  Square  Law 

dx  =  -a(t)-y(t)  dt;  (7) 

dy  =  -b(t) • x(t)  dt. 
Assume  f  is  1-1  and  dif f erentiable  so  that  dT  =  f'(t)dt  and  t  = 
f_1(T).   Substituting  for  t  and  dt  in  equations  (7)  gives 

dCxff^CT))]  =  -a(f  1(T))-y(f  1(T))-dT/f ■  (t)  (8) 

=  dx(f  1(T))/f '  (t) 
(assuming  f'(t)  *  0)  and 


d[y(f'1(T))]  =  -b(f1(T))-x(f  1(T))-dT/f  ■  (t)         (9) 

■-1 


=  dy(f-'(T))/f  (t). 


Letting 


X(T)  =  x(f1(T))  , 
Y(T)  =  y(f1(T)), 
A(T)  =  a(f  1(T))/f  '  (f'1(T)), 
B(T)  =  b(f1(T))/f'(f1(T)), 
it  follows  that 

dX(T)  =  -A(T)-Y(T)  dT;  (10) 
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dY(T)  =  -B(T) • X(T)  dT, 
which  is  just  the  Lanchester  Square  Law  with  a  new  time  parameter 
T  and  "transformed"  attrition  coefficients. 
Now  suppose 

c 

fit)    =  a.'jb{*t)dx  (11) 

0 

so  f'(t)  =  ab(t).   Then  f'(f"1(T)  =  ab(f"1(T),  so  B(T)  =  a.   With 
the  transformation  defined  in  equation  (11)  the  attrition 
coefficient  B  in  equations  (10)  is  made  constant.   This  implies 
one  can  assume  without  loss  of  qenerality  that  b  is  constant  in 
equations  (7) . 
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Appendix  2.  Variance  Reduction 

We  consider  here  a  method  of  computing  Ay  and  Ax  within  a 
time  interval  At(  which  provides  a  reduction  in  the  variance  of 
these  statistics.   In  the  Lanchester  Square  Law,  these  quantities 
represent  the  numbers  of  casualties  occurring  on  sides  Y  and  X, 
respectively,  during  the  time  interval  in  question.   But  within 
the  simulated  battles  (whether  at  the  NTC  or  with  the  Janus 
combat  simulation) ,  the  numbers  of  casualties  are  outcomes 
depending  on  numerous  random  events  and  conditions  occurring  as 
the  battle  evolves  in  time.   We  thus  consider  determination  of  Ay 
and  Ax  as  a  statistical  estimation  problem.   We  suppose  there  are 
true  (but  unknown)  mean  values  of  the  distributions  of  numbers  of 
casualties  on  each  side  within  the  ith  time  interval,  and  we 
consider  the  task  of  estimating  these  means.   The  estimates  will 
then  be  used  as  input  values  of  Ay  and  Ax  in  the  computation  of 
the  battle  trace. 

An  obvious  way  to  estimate  the  expected  number  of  battle 
units  of  a  certain  type  that  would  be  killed  in  a  future  battle 
with  identical  conditions  is  to  count  the  number  of  units  of  this 
type  that  were  killed  in  the  simulated  battle  that  is  under 
consideration.   This  estimator,  known  as  the  body  count 
estimator,  is  not  commonly  used  in  practice  because  a  superior 
estimator,  the  "sum  of  pk's"  estimator,  is  available.   Suppose  M 
engagements  occur  in  the  ith  time  interval,  and  the  kth  of  these 
has  kill  probability  pk.   A  sum  of  pk's  estimate  is  formed  by 
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summing  the  pk '  s  used  in  the  simulated  battle  over  all 
engagements  against  targets  of  the  type  in  guestion. 

The  body  count  and  sum  of  pk's  estimators  can  be  compared  as 
follows.   Suppose  results  of  engagements  against  targets  of  the 
type  in  question  are  represented  by  Bernoulli  random  variables 
X1 ,  X2,  ...,  XM.   Here,  Xk  is  zero  if  the  kth  engagement  against  a 
target  of  the  given  type  does  not  result  in  a  kill  (which  will 
happen  with  probability  1  -  pk,  say) ,  and  Xk  is  1  if  the  kth 
engagement  does  result  in  a  kill.   The  random  variable  M  is  the 
number  of  engagements  against  the  type  of  target  in  question  in 
the  simulated  battle.   The  number  of  casualties  of  the  given  type 
is  Sk=1MXk  so  the  body  count  estimator  is  C  =  Sk=1MXk.   The  sum  of 
pk's  estimator  is  C*  =  Zk=1Mpk.   Since  the  expected  value  of  Xk  is 
pk,  both  C  and  C*  have  expected  value  E[Sk=1Mpk]  ,  where  "E"  denotes 
expectation  with  respect  to  the  distribution  of  M.   This  is  the 
correct  mean  (that  is,  both  C  and  C*  are  unbiased  estimators). 
However,  the  variance  of  C  is 

E[2k=1Mpk(l-pk)]  +  V[Sn=1M  pk], 
where  "V"  denotes  variance  with  respect  to  the  distribution  of  M. 
The  variance  of  C*  is  equal  to  only  the  second  term  in  this 
expression  and  is  thus  smaller  than  the  variance  of  C.   In 
statistical  terminology,  the  sum  of  pk  •  s  estimator  is  better  than 
the  body  count  estimator.   As  mentioned  above,  the  kth  term  in 
the  sum,  pk,  is  the  expected  number  of  kills  on  the  kth 
engagement,  since  the  expected  value  of  Xk  is  pk.   Thus  the  sum 
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of  pk ' s  estimator  can  be  interpreted  as  summing  fractional 
expected  kills  resulting  from  each  engagement. 

Appendix  3.  Approximate  Variance  of  BL 

It  may  be  desirable  to  plot  an  indicator  of  the  standard 
error  of  the  battle  trace  along  with  the  trace  values  themselves. 
In  what  follows  we  consider  how  the  variance  of  the  log  trace 
might  be  approximated. 

As  indicated  in  Section  1,  the  log  trace  is  defined  by 
BL=[ln(Ay+a)  +  ln(x+a)]  -  [ln(Ax+a)  +  ln(y+a)],        (12) 
where  a  is  a  suitable  positive  constant.   The  specific  value  of  a 
has  relatively  little  effect  on  our  considerations  here,  so  we 
ignore  it  in  the  remainder  of  this  Appendix.   Assume  the 
quantities  associated  with  force  sizes  in  the  second  Lanchester 
Square  Law  equations  (7)  are  random;  these  equations  suggest  -Ay 
and  x  have  positive  covariance,  and  similarly  for  -Ax  and  y. 

Now  let  us  consider  the  variance  of  [ln(Ay)  +  ln(x)].    Take 
the  log  of  both  sides  of  the  second  equation  in  (7)  to  obtain 

ln(-Ay)  =  ln(b)  +  ln(x)  +  ln(At), 
or,  assigning  simplifying  symbols, 

Y  =  U  +  X  +  k, 
where  k  is  constant  (or  at  worst  is  not  correlated  with  X) . 
Assume  temporarily  that  the  random  variables  have  been  centered 
at  their  means,  so 

Cov(Y,X)  =  E(Y- X) 
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=  E( [U  +  X  +  k]X) 

=  E(UX  +  X2  +  kX) 

=  V(X), 
assuming  U  and  X  are  uncorrelated. 
It  follows  that 

V(Y  +  X)  =  V(Y)  +  V(X)  +  2Cov(Y,X) 

=  V(Y)  +  3V(X) 

=  V(U)  +  V(X)  +  3V(X) 

=  V(U)  +  4V(X) . 
With  a  similar  argument  for  the  second  term  in  expression  (7), 
the  variance  of  BL  is,  under  our  assumptions, 

V(BL)  =  4[V(ln(x))  +  V(ln(y))]  +  [V(ln(b))  +  V(ln(a))].    (13) 
Since  we  anticipate  that  ln(b)  and  ln(a)  will  have  small 
variances  relative  to  those  of  ln(x)  and  ln(y),  it  may  reasonable 
to  approximate  V(BL)  by  the  first  term  in  expression  (13). 

The  variances  of  ln(x)  and  ln(y)  can  be  estimated  through  a 
sample  of  their  respective  values  in  the  preceding  few  time 
intervals  (as  determined  by  the  sum  of  maxp  estimators) .   As  an 
alternative,  since  we  plan  to  use  the  estimated  variances  to 
establish  a  measure  of  the  guality  of  the  estimated  battle  trace, 
it  may  suffice  to  use  a  quantity  that  only  "tracks"  the 
underlying  standard  deviation.   The  variances  of  ln(x)  and  ln(y) 
can  be  roughly  bounded  as  follows.   Following  the  argument  for 
the  variance  of  the  sum  of  pk '  s  estimator  C*  above,  it  can  be 
seen  that  the  sum  over  X-units  of  maxpx  (see  Appendix  2)  has 
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variance  approximately  equal  to 

EN[£x-unitsmaxPx  I  Nl  ~  N-maxpx,  (14) 

where  "EN"  denotes  the  expectation  with  respect  to  the  random 
number  N  of  X-units  involved  in  the  ith  time  interval.   Since  x 
will  tend  to  be  greater  than  1,  ln(x)  will  have  variance  less 
than  that  of  x.   Thus  a  bounding  value  may  be  obtained  through 
expression  (14),  as  follows.   The  averages  on  the  right-hand  side 
of  (14)  can  either  be  estimated  within  the  current  time  interval 
by  using  the  observed  values  n  and  maxpx,  or  a  pooled  estimate 
can  be  maintained,  combining  with  data  from  previous  time 
intervals.   It  might  be  noted  that  the  first  of  these  estimators 
suggests  that,  very  roughly,  the  standard  deviation  of  BL  may 
parallel  (fluctuate  with)  2/(x  +  y) ,  or  for  our  purposes,  simply 
/(x  +  y) .   Some  experimentation  will  be  necessary  to  see  whether 
this  value  gives  a  useful  indication  of  the  variation  seen  in  BL . 

Appendix  4.   Plots  of  Log  Trace  for  Janus  Runs 

Example  log  traces  were  made  with  Janus  runs  using  the 
implementation  described  in  Section  5  and  documented  in  Appendix 
5.   The  runs  all  used  a  "Fulda  Gap"  scenario  developed  at  TRAC- 
Monterey  as  part  of  a  Janus  Tutorial.   In  this  scenario  a  Blue 
force  defends  against  a  numerically  superior  Red  attacker.   The 
tutorial  concerns  placement  of  a  Blue  task  force  (which  we  call 
"task  force  Blue  1")  consisting  primarily  of  tanks  and  armored 
personnel  carriers  overlooking  a  river  crossing  that  is  likely  to 
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be  used  by  the  attacker.   We  consider  also  a  part  of  the  Red 
force  that  eventually  engages  with  Blue  1  to  be  a  Red  task  force 
("Red  1") .   We  ran  seven  iterations  of  this  scenario;  several  of 
the  battle  trace  plots  that  follow  show  plots  for  all  seven  runs 
superimposed  on  one  figure.   The  intent  of  these  displays  is  to 
illustrate  how  battle  traces  can  vary  from  one  iteration  of  a 
battle  to  another. 

Most  of  the  figures  relate  to  the  Blue  1  task  force  or  to 
the  Red  force.   In  some  cases  the  Red  1  task  force  is  considered. 
These  figures  suggest  Blue  1  enjoyed  success  early  in  each 
battle,  but  for  times  beyond  about  two  minutes  the  Red  force  is 
generally  winning.   In  several  of  the  figures  it  can  be  seen  that 
Blue  1  suffers  particular  difficulties  around  time  t=4  minutes. 
This  corresponds  to  a  point  in  the  scenario  at  which  surviving 
units  of  Blue  task  force  1  are  set  in  motion  for  the  first  time 
in  the  battle.   The  relevant  battle  traces  suggest  this  movement 
is  initiated  too  far  into  the  battle,  and  it  proves  disastrous 
for  Blue  1. 
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Battle  Trace,  Fulda  Gap  Scenario 

Blue  Side,  Smoothed  (4  pt) 
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Figure  1.   Smoothed  log  trace  for  Blue  1  on  one  iteration  of  the 
Fulda  Gap  scenario.   It  can  be  seen  that  Blue  1  does  well 
initially  (BL>0)  but  beyond  t=2  minutes  of  battle,  the  Red  forces 
are  generally  winning.   Tick  marks  along  the  bottom  of  the  figure 
correspond  to  times  at  which  points  on  BL  were  computed.   The 
threshold  used  is  32  (i.e.,  a  point  is  calculated  on  BL  when 
approximately  32  shots  have  been  fired  by  or  at  Blue  1. 
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Battle  Trace,  Fulda  Gap  Scenario 


Red  Side 


3  4 
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Figure  2.   Non-smoothed  log  trace  for  all  Red  forces, 
corresponding  to  the  iteration  shown  in  Figure  1  (threshold  = 
32).   Note  the  general  improvement  in  Red's  log  trace  over  time; 
note  also  the  variations  in  BL  from  one  time  interval  to  the 
next. 
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Figure  3.   Blue  1  log  traces  for  seven  iterations  of  the  Fulda 
Gap  scenario  (threshold  =  32) .   Note  strong  similarity  of  the  log 
traces  in  overall  shape,  despite  variations  between  iterations. 
The  general  sharp  drop  in  BL  near  t=4  minutes  corresponds  to  a 
point  in  the  scenario  at  which  Blue  1  starts  to  move  its 
surviving  units.   Units  on  axes  appear  fuzzy  because  multiple 
overlays  were  used  to  construct  the  figure. 
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Figure  4 .   Red  task  force  1  log  traces  corresponding  to  the  seven 
iterations  shown  in  Figure  3. 
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Battle  Trace,  Fulda  Gap  Scenario 


Time  (min) 


-^*—  Log  Trace 


CER 


Figure  5.   Log  trace  for  Blue  1,  together  with  plot  of  casualty 
exchange  ratios  (CERs)  in  each  time  interval.   Note  similarity  of 
shapes  of  the  two  plots. 
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log  BT  vs  log  CER 

threshold  32 


CT> 
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Observed  Pairs 


Line  with  Slope  2 


Figure  6.   Log-log  plot  of  battle  trace  versus  CER.   Note  line 
with  slope  2  provides  good  fit  of  the  plotted  points. 
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Attrition  Ratio  vs.  Time 

threshold  32 


Time  (min) 


Figure  7.   Plot  of  Blue  1  versus  Red  force  attrition  coefficient 
ratio  (B/a)  estimates.   The  values  were  estimated  through 
residuals  relative  to  the  line  in  Figure  6. 
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Figure  8.   Log  Blue  1  versus  Red  force  attrition  coefficient 
ratio  estimates,  ln(B/a),  plotted  as  a  function  of  time  into  the 
battle  for  seven  iterations  of  the  Fulda  Gap  scenario.   These 
iterations  correspond  to  the  battle  traces  shown  in  Figure  3. 
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Figure  9.   Plots  of  Bc  for  seven  iterations  of  the  Fulda  Gap 
scenario  corresponding  to  the  battle  traces  shown  in  Figure  3. 
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Appendix  5.   Fortran  Code  for  Battle  Trace  in  Janus 

Implementation  of  the  battle  trace  in  Janus (A)  at  TRAC- 
Monterey  required  making  additions  to  several  subroutines  within 
the  Janus  Fortran  source  code.   We  include  listings  of  the  major 
subroutines  affected  for  the  purpose  of  documentation.   As  may  be 
seen,  the  additions  have  been  commented  rather  generously,  so 
they  are  "self  documented"  to  a  large  extent. 
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C SUBROUTINE- -RELOAD 

SUBROUTINE   RELOAD 


A.D.KELLNER,  TRAC-WSMR 


C*****   MODIFIED  FOR  BATTLE  TRACE  COMPUTATION  *****  J.C.  HOFFMAN,  TRAC-MTRY 
C  D.R.  BARR,  NPS 

C  2  APR  92 


C- 

c 

c 

c 

c 

c 

c 

c 

c- 

c 

c 

c 

c 

c 

c 

c- 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c- 


PURPOSE: 


INPUTS: 


Sub-driver  to  process  a  Direct  Fire  event.  This  subrou- 
tine performs  TARGET  SELECTION  for  a  unit  firing  a 
Direct  Fire  weapon;  calls  "SHOOT"  to  simulate/evaluate 
the  firing  event;  then  calls  "RLNEXT"  to  schedule  the 
next  Direct  Fire  event  to  be  processed. 


(  From  COMMON  /LOAD/   in   GLOBLOAD . FOR  ) 


NXTUNIT   -   Unit  number  of  shooter 
NXTSIDE   -   Side  of  shooter 


DICTIONARY: 


IUNIT,  IS 
JUNIT,  JSIDE 

IWPN 
IWPNSLT 


Unit  number,  side  of  shooter 
Unit  number,  side  of  target 

Absolute  Weapon  type  (1-1 
System  Relative  Weapon  (1-10) 


INCLUDE 
INCLUDE 
INCLUDE 
INCLUDE 
INCLUDE 
INCLUDE 
INCLUDE 
INCLUDE 

INCLUDE 


JGLOBE: GLOBAL. F 
JGLOBE : GLOBUNITS . FOR 
JGLOBE : GLBMOPPS . FOR ' 
JGLOBE : GLBWEAPN . FOR ' 
JGLOBE : GLOBFI R . FOR ' 
JGLOBE : GLOBACQ . FOR ' 
JGLOBE : GLOBLOAD . FOR ' 
JGLOBE : GLOBMOVE . FOR ' 

JGLOBE : GLBTRACE . FOR ' 


!  CLOCK 

!  KSYSTYP, FIREPS, KFIR1 

!  KGUIDE 


I  IFLYMODE 

!  BATTLE  TRACE  Variables*** 


DIMENSION 
DIMENSION 


IENEM(NMVISB) ,  PKS(NMVISB),  RNGTAR (NMVISB) 
IWPNSAV(NMVISB) ,  ISTATS (NMVISB ) ,  VALUE (NMVI SB ! 


IPRINT  =  0 


IUNIT 
ISIDE 


NXTUNIT 
NXTSIDE 


Make  popped  special-special  flyers  go  through  SFRELOAD  

IF(  ISIDE  .EQ.  1  .AND. 

FLYERS (KSYSTYP( IUNIT, ISIDE) , ISIDE)  . EQ .  NUMFLYRS   .AND 
IFLYMODE (IUNIT, ISIDE)  . LT .  0  )  THEN 

CALL  SFRELOADf IUNIT, ISIDE) 
GOTO  999 
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ENDIF 


IAMFIRNGdUNIT,  ISIDE)   =   0 
FIRNXT (IUNIT, ISIDE)   =   99999 


Clear  unit's  "Firing"  flag 


D 
D 
D 
D 
D 
D 
D 
D 

D 
D 
D 

D 
D 
D 
D 
D 
D 
D 
D 
D 
D 
D 
D 


IF(  IFIRST  .NE.  99  )   THEN 

IFIRST   =   99 

TYPE  * 

TYPE  *, 'Enter  RELOAD  unit  number  for  debug  print:' 

ACCEPT  *,  I  PUN IT 

TYPE  *, 'Enter  RELOAD  side  for  debug  print:' 

ACCEPT  *,  IPSIDE 

ENDIF 

I  PR  INT  =  0 

IF(  ( IUNIT. EQ. I  PUN IT)  .AND.  ( ISIDE . EQ . IPSIDE)  )   IPRINT  =  99 

IF(  (IPUNIT  .EQ.  0  )   .AND.  ( ISIDE . EQ . IPSIDE)  )   IPRINT  =  99 


IF(  IPR 
TYPE  * 
TYPE  *, 
TYPE  *, 
TYPE  *, 
TYPE  * 
ITYPE 
TYPE  *, 
TYPE  * 
TYPE  *, 
TYPE  *, 
END  IF 


INT  .GT.  0 


THEN 


SUBROUTINE  RELOAD 


KSYSTYP) IUNIT, ISIDE) 
■-  GAME  TIME  (Sec)  =',  CLOCK*60.0 

■-  IUNIT, ISIDE, ITYPE  =',  IUNIT, ISIDE, ITYPE 

•-  IAMFIRNGdUNIT,  ISIDE)  =',  IAMFIRNG  (  IUNIT  ,  ISIDE) 


IF(  IDEFLdUNIT,  ISIDE)   .  EQ .   2 

IF(  IDEFLdUNIT,  ISIDE)   .  EQ  .  -2 

IF(  NSCORE( IUNIT, ISIDE)  . LT .  1  ) 

IF(  KINOPSTATt IUNIT, ISIDE)  . GT . 
ITYPE   =   KSYSTYP (IUNIT, ISIDE) 

IF(  FIRERS( ITYPE, ISIDE)  . EQ .  2  ) 

IF(  FIRERS( ITYPE, ISIDE)  . LE .  0  ) 


GOTO  3  i0 
GOTO  300 
GOTO  300 
GOTO  300 

GOTO  3  00 
GOTO  3  00 


IF(  MOUNTED ( IUNIT, ISIDE)  .GT.  0  )   THEN 


Unit  is  a  passenger,  update  his  current  status 
IHOST   =   MOUNTEDdUNIT,  ISIDE) 


XUNIT( IUNIT, ISIDE) 
YUNIT( IUNIT, ISIDE) 
SPDU( IUNIT, ISIDE) 
IHOLFIR( IUNIT, ISIDE)  = 
TSUPRSt IUNIT, ISIDE) 


XUNIT( IHOST, ISIDE) 
YUNITdHOST,  ISIDE) 
SPDU( IHOST, ISIDE) 
IHOLFIR( IHOST, ISIDE] 
TSUPRS( IHOST, ISIDE) 


Full  defilade 
Full  defilade 
Dead 

CHEM 


ENDIF 


IF(  IHOLFIRdUNIT,  ISIDE; 

FIRNXT ( IUNIT, ISIDE) 

GOTO  3  00 
ENDIF 


NE.  0  )  THEN 
CLOCK  +0.2 
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IF(  TSUPRSf IUNIT, ISIDE)  .GT.  CLOCK  )   THEN 

FIRNXTdUNIT,  ISIDE)   =   TSUPRS  (  IUNIT,  ISIDE) 

GOTO  3  00 
END  IF 

FIRNXTdUNIT,  ISIDE)   =   99999. 
JSIDE   =   3  -  ISIDE 

JUNIT    =   0 
SUMVAL   =   0 . 
ITGT     =   0 

ITASK    =   KTASKFOR( IUNIT, ISIDE)  !  B-Trace  MOE  *** 

C  ITask  Force  of  Observer 

D  IF(  I  PR  I NT  .GT.  0  )   THEN 

D  TYPE  * 

D  TYPE  * ,  ' LOOK  FOR  A  TARGET 

D  ENDIF 

DO  100   NM  =  1,  NM7ISB 

IF(  KDETECTD(NM, IUNIT, ISIDE)  . NE .  1  )   GOTO  100      !  Not  detected 

JUNIT   =   NMYUNIKNM,  IUNIT,  ISIDE) 

IF(  JUNIT  .EQ.  0  )   GOTO  I  !  Empty  slot 

JTYPE     =   KSYSTYP (JUNIT, JSIDE) 
JFIRPRI   =   KFIRPRI (ITYPE, JTYPE, ISIDE) 

JTASK     =   KTASKFOR {JUNIT, JSIDE)  IB-Trace  MOE  *** 

ITask  Force  cf  the  target 

IF(  JFIRPRI  .LE.  0  )   GOTO  100  !  No  priority 

IF(  NSCORE (JUNIT, JSIDE)  . LT .   1  )   GOTO  100  !  Already  dead 

FIRNXTdUNIT,  ISIDE)   =   CLOCK  +0.1 

C Fetch  Target  range 

DX  -      XUNTT (IUNIT, ISIDE)  -  XUNIT ( JUNIT , JSIDE)        !  From  target 
DY   =   YUNITI IUNIT, ISIDE)  -  YUNIT (JUNIT, JSIDE)        !  to  firer 
RANGE   =   SORT(  DX*DX  +  DY*DY  ) 

C Check  if  Target  beyond  max  firing  range 

RNGMAX   =   PRIMRNG( ITYPE, ISIDE) 
IF(  RANGE  .GT.  RNGMAX  )   GOTO  100 


C Check  if  Target  in  full  defilade 

IF(  IDEFL(JUNIT, JSIDE)   . EQ .   2   .OR. 
*  IDEFL(JUNIT, JSIDE)  . EQ .  -2         )   THEN 

C Don't  fire  unless  within  50  meters 

IF(  RANGE  .GT.  0.050  )   GOTO  100 

ENDIF 
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IF(  I  PRINT  .GT.  0  )   THEN 


D 

TYPE  * 

D 

TYPE  * 

D 

TYPE  * 

D 

TYPE  * 

D 

TYPE  * 

D 

END  IF 

POSSIBLE  TARGET:' 

■--  JUNIT,  JSIDE   =',  JUNIT, JSIDE 
■--  JTYPE,  JPRI    =',  JTYPE, JFIRPRI 
■--  RANGE,  RNGMAX  =',  RANGE,  RNGMAX 


C Fetch  weapon  to  use  against  this  target 

C Force  special -special  flyers  not  to  use  rf  missile  by  setting 

C  rounds  left  to  zero 

IFLYTYPE   =   FLYERS ( ITYPE, ISIDE) 

IF(  ISIDE. EQ.l  .AND.  IFLYTYPE . EQ . NUMFLYRS  )   THEN 

ISAVRNDS   =   KRLEFTU,  IUNIT,  ISIDE) 

KRLEFTd,  IUNIT,  ISIDE)   =   0 
END  IF 


CALL  WTCHWPN  (  IUNIT, ITYPE, ISIDE, JTYPE, RANGE,  IWPNSLT, IWPN  ) 

D  IF{  I  PRINT  .GT.  0  )   THEN 

D  TYPE  *, ' IWPNSLT, IWPN  =',  IWPNSLT, IWPN 

D  END  IF 


Restore  spec ia 1 -spec ia 1  rf  missile  rounds  previously  zeroed. 


IF(  ISIDE. EQ.l  .AND.  IFLYTYPE . EQ . NUMFLYRS  )   THEN 

KRLEFTd,  IUNIT,  ISIDE)   =   ISAVRNDS 
END  IF 


Break  if  no  weapon  selected 
IF(  IWPNSLT  .LE.  0  )   GOTO  100 


c Get  SSKP  for  this  target 

CALL  PROB  (  IUNIT,  IWPN,  JUNIT,  ISIDE,   JSIDE, 

RANGE,   DY,    DX,    ISTAT,   SSKP    ) 

IF{  MOPPdUNIT,  ISIDE)  .  EQ .  1  ) 
*  SSKP   =   SSKP  *  PHMOPPdWPN,  ISIDE)        !  CHEM 

D  IF(  I  PRINT  .GT.  0  )   THEN 

D  TYPE  *, ' SSKP  =',  SSKP 

D  END  IF 

IF(  SSKP  .LT.  0.0  5  )   THEN 

IF(  SSKP  .GT.  0.0  )   FIRNXTt IUNIT, ISIDE)  =  CLOCK  +0.02 

GOTO  100 
ENDIF 
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q Check  for  guided  weapon  that  can't  track  thru  smoke 

IF(  KGUIDEdWPN,  ISIDE)  .  EQ .  2  )  THEN 

CALL  UNITXYZ  (  IUNIT, ISIDE, ITYPE,  XF,YF,ZF  ) 

JTYPE   =   KSYSTYPfJUNIT, JSIDE) 

CALL  UNITXYZ  (  JUNIT, JSIDE, JTYPE,  XT,YT,ZT  ) 

ISENS   =   1  !  OPTICAL 

DX   =   ABS(  XF-XT  ) 
DY   =   ABS(  YF-YT  ) 

IF(  (DX.GT. 0.025)  .OR.  ( DY . GT . 0 . 02 5 )  )   THEN 

CALL  SMOKELOS  (  XF,YF,ZF,  XT,YT,ZT,  ISENS,  ISEE  ) 
IF(  ISEE  .NE.  1  )  GOTO  100 

CALL  DOLASLOS  (  XF,YF,ZF,  XT,YT,ZT,  ISENS,  OLEN  ) 
IF(  OLEN  .GT.  2.0  )  GOTO  100 

END  IF 

ENDIF 

C Add  potential  target  to  local  target  list 

ITGT  -       ITGT  +  1 

IENEM(ITGT)  =  JUN 

RNGTAR(ITGT)  =  RANGE 

IWPNSAV(ITGT)  =  IWPNSLT 

ISTATS(ITGT)  =  ISTAT 

PKS(ITGT)  =  SSKP 

VALUE (ITGT)  =  SSKP  *  FLOAT (  JFIRPRI  ) 

SUMVAL  =  SUMVAL  +  VALUE (ITGT) 

c*****    UPDATE  BATTLE  TRACE  1BTRACE 

C  See  if  there  is  a  new  maximum  sspk  value  for  each        !BTi 

C  enemy  task  force.   If  so,  then  record  the  value  !BT; 

C  in  the  array  BTRACEVAL  1BTRACE 

IF  (  SSKP  .GT.  BTRACEVAL {ISIDE, IUNIT, JTASK)  )  THEN  ! BTRACE 

BTRACEVAL (ISIDE, IUNIT, JTASK)  =  SSKP         ! BTRACE 

ENDIF  ! BTRACE 

c****«    END  0F  BATTLE  TRACE  UPDATE  i BTRACE 

100    CONTINUE 


Skip  if  no  valid  targets 
IF(  ITGT  .EQ.  0  )   GOTO  300 

If  only  one  target,   select  it 

IF(  ITGT  .EQ.  1  )   THEN 

JUNIT     =  IENEM(ITGT) 

SSKP      =  PKS(ITGT) 

RANGE     =  RNGTAR(ITGT) 

IWPNSLT   =  IWPNSAV(ITGT) 
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I  STAT 
GOTO  2  00 
END  IF 


istatsqtgt: 


Target  with  highest  value  (  SSKP  *  JFIRPRI 
has  higest  probability  of  being  selected 


CUMRAT  =0.0 

CALL  UN I RAN  (DRAW) 

CHECK  =  DRAW  *  SUMVAL 

IF(  I  PRINT  .GT.  0  )   THEN 

TYPE  * 

TYPE  *, 'DRAW,  CHECK  =',  DRAW, CHECK 

END  IF 


DO  150   I  : 
CUM PAT 


1,  ITGT 

=   CUMRAT  +  VALUE { I ) 


D 
D 

D 


150 


IF(  I  PRINT  -GT.  0  )   THEN 

TYPE  *, ' LOOP  I,  CUMRAT  =' 

END  IF 


IF(  CUMRAT  .LE.  CHECK  ) 

JUNIT    =  IENEM(I) 

SSKP     =  PKS(I) 

RANGE    =  RNGTAR(I) 

IWPNSLT  =  IWPNSAV(I) 

ISTAT    =  ISTATS(I) 
GOTO  2  00 


I , CUMPAI 
GOTO  15  0 


NTINUE 


200 


A  target  has  been  selected 
CONTINUE 


D 
D 
D 
D 
D 


*  ***** 


IF(  IFF  INI  .GT.  0  )   THEN 

TYPE  * 

TYPE  *, 'RELOAD  HAS  SELECTED  TARGET  UNIT  =  ',  JUNIT 

TYPE  *,' SSKP  =',  SSKP 

END  IF 

UPDATE  BATTLE  TRACE  


Update  the  influence  of  all  shots  fired  on  both 
friendly  and  enemy  taskforces. 


iBTRACE 

IBTRACE 
IBTRACE 


Influence  of  the  shot  on  the  TARGET  Task  Force  IBTRACE 

SIPCTPKS(JSIDE, JTASK)  =  SIPCTPKS ( JSIDE, JTASK)  +  SSKP  IBTRACE 

KNUMSHOTS (ISIDE, ITASK)  =  KNUMSHOTS ( ISIDE, ITASK)  +  1  IBTRACE 

KNUMSHOTS (JSIDE, JTASK)  =  KNUMSHOTS (JSIDE,  JTASK )  +  1  IBTRACE 

The  array  KNUMSHOTS  is  zeroed  each  time  the  Battle  IBTRACE 

is  computed  for  a  given  task  force.   This  is  done  IBTRACE 

in  subprogram  BTCOMP  called  in  RUNJANUS  IBTRACE 
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2  1 1) 


Influence  of  the  shot  on  the  FIRER's  Total  Force  ! BTRACE 

DO  210  I  =  1,  NUMTASKS  ! BTRACE 

IF  (  BTINFL(JSIDE, JUNIT, I)   .GT.  0  )  THEN  ! BTRACE 

SSHOTPKSdSIDE,  I)  =  SSHOTPKS  (  ISIDE,  I  )  +  SSKP  !  Bl  I 

END IF  ! BTRACE 

CONTINUE  !  BTRACE 


c*****   END  of  BATTLE  TRACE  UPDATE 
IAMFIRNGdUNIT,  ISIDE)   =   1 


! BTRACE 
!  Set  unit's  "Firing"  flag 


CALL  SHOOT  (  IUNI1,  ISIDE,  ITYPE,  IWPNSLT, 

JUNIT,  SSKP,   RANGE,  ISTAT,  KILLS   ) 


3  00    CONTINUE 

Schedule  next  firing  event  &  RETURN 

CALL  RLNEXT 


IF(  IPRIN1 
TYPE  * 
TYPE  *,'-- 
TYPE  *, 

TYPE  *, 
TYPE  *, 
TYPE  * 
TYPE  * 
END  IF 


.GT.  0  ) 


THEN 


EHD  RELOAD  

-  FIRNXT    =',  FIRNXTdUNIT,  ISIDE)  * 

-  IAMFIRNG  =',  IAMFIRNGdUNIT,  ISIDE) 

-  TNEXT(4)  =',  TNEXT(4)  *  60. 


999 


CONTINUE 


RETURN 
END 
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C SUBROUTINE- -BTCOMP- 

C 
C 

SUBROUTINE  BTCOMP 


-D.R.  BARR,  DEPT  OF  MATH,  MPS 
J.C.  HOFFMAN,  TRAC-MTRY 
6  APR  92 


q This  subroutine  sets  up  Battle  Trace  data  structures  needed 

C        to  compute  the  Battle  Trace  for  all  Janus (A)  Task  Forces. 


C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C 
C* 


********************* 


LOCAL  VARIABLES 


********************************* 


ISIDE  - 
IUNIT  - 
ITASK  - 
JTASK  - 
PKTASKMAX  - 

PKUNITMAX  - 


Index  of  friendly  side 

Index  of  friendly  unit 

Index  of  friendly  task  force 

Index  of  enemy  task  force 

The  maximum  single  shot  kill  probability 
recorded  against  each  enemy  task  force 

The  sum  of  BTRACE  over  all  unit,  a:   .gned  to 
a  given  friendly  task  force 


***************************** 


****** 


*********** 


*  *  *■  -x   * 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 
***c 
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50 


INCLUDE 
INCLUDE 
INCLUDE 

REAL 


' JGLOBE : GLOBAL . FOR ' 

' JGLOBE : GLOBUNITS . FOR ' 

' JGLOBE : GLBTRACE . FOR ' 


Battle  Trace  Variables 


PKTASKMAX ( NUMS I DES,NUMT ASKS) , 
PKUNITMAX (NUMS IDES, NUMT ASKS ) 


Zero  influence  accumulators 

DO  50  I  =  1,  NUMS IDES 
DO  4  0  J  =  1,  NUMT ASKS 
PKTASKMAX (I, J)  =0.0 
PKUNITMAX (I, J)  =0.0 
CONTINUE 
CONTINUE 


Examine  each  value  for  the  maximum  SSKP  recorded  for  each 
unit  on  each  side.   Determine  the  maximum  SSKP  recorded  for  all 
enemy  task  forces.   This  value  represents  the  contribution  of 
each  particular  unit  to  the  battle. 

Next,  compute  the  task  force  contributions  for  both  friendly  and 
enemy  task  forces . 

a.  Sum  the  contributions  for  each  friendly  task  force 
which  has  reach  the  THRESHOLD  for  Battle  Trace  computation  and 
reinitialize  the  BTMAX  values  to  zero  (This  last  action  has  the 
effect  of  starting  a  new  computation  cycle  for  the  Battle 
Trace. ) 

b.  Compute  the  contribution  of  all  friendly  units  which 
detect  at  least  one  element  of  each  enemy  task  force  for  those 
enemy  task  forces  which  have  reached  the  THRESHOLD  for  battle  trace 
computation.   Do  this  by  summing  BTRACEVAL  for  all  friendly  units 
indexed  by  the  appropriate  enemy  task  force.   These  sums  are  stored 
in  array  PKUNITMAX.   Reinitialize  the  appropriate  entries  in 
BTRACVAL  by  setting  to  zero  all  entries  which  correspond  to  units 
which  contribute  to  the  computation  of  battle  trace  for  a 
particular  enemy  task  force. 
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DO  3  00  ISIDE  =  1,  NUMSIDES 

JSIDE  =  3  -  ISIDE 

DO  200  IUNIT  =  1,  KNUMUNITSt ISIDE) 

DO  100  JTASK  =  1,  NUMTASKS 

C Determine  maximum  SSKP  for  each  friendly  unit  against 

C        all  enemy  task  forces 

IF  (  BTRACEVAM ISIDE, IUNIT, JTASK)  .GT. 

*  BTMAX( ISIDE, IUNIT)  )  THEN 
BTMAX( ISIDE, IUNIT)  =  BTRACEVAL ( ISIDE, IUNIT, JTASK ) 

END  IF 

C Compute  the  contribution  for  friendly  units  which  observe  at  least 

C  one  unit  of  an  enemy  task  force  for  those  enemy  task  forces  which 
C        have  reached  the  threshold  for  battle  trace  computation. 

IF  (  KNUMSHOTS (JSIDE, JTASK)  .GE.  THRESHOLD  )  THEN 

PKUNITMAX (JSIDE, JTASK)  =  PKUNITMAX (JSIDE, JTASK)  + 

*  BTRACEVAL (IS IDE, IUNIT, JTASK) 

C Reinitialize  the  contribution  for  each  friendly  unit  which 

C        contributed  to  the  calculation  of  battle  trace  for  an  enemy 
C        task  force. 

BTF     \L( ISIDE, IUNIT, JTASK)  =  0 
ENDIF 

100        CONTINUE 

ITASK  =  KTASKFOR( IUNIT, ISIDE) 

C Sum  the  maximum  contributions  to  the  battle  trace  if  the  battle 

C        trace  is  to  be  computed  for  a  friendly  task  force. 

C C 

C  Compute  the  max  reduction  over  columns  of  BTRACEVAL.  PKTASKMAX  C 
C  now  contains  the  maximum  single  shot  kill  probability  of  all  C 
C  units  of  a  particular  friendly  task  force  who  detect  (and  C 
C  have  a  possibility  of  firing  upon)  at  least  one  unit  of  an  C 
C        opposing  task  force  C 

c c 

IF  (  KNUMSHOTS (ISIDE, ITASK)   . GE .  THRESHOLD  )  THEN 
PKTASKMAX (ISIDE, ITASK)  =  PKTASKMAX ( ISIDE, ITASK)  + 

*  BTMAX( ISIDE, IUNIT) 

C Reinitialize  max  contribution  for  next  computation  interval 

C  (NOTE:   The  interval  of  computation  in  terms  of  battle  time 

C  is  a  function  of  the  dynamics  of  the  scenario  and  the  value 

C  set  in  THRESHOLD.   Larger  values  of  THRESHOLD  will  result 

C  in  a  longer  period  of  battle  time  between  subsequent  calculation 

C  of  the  battle  trace.   Likewise,  battle  periods  which  do  not  contain 

C  a  sufficient  amount  of  weapon  firings  or  shots  received  by  a 

C  particular  task  force  will  not  attain  the  THRESHOLD  for  battle 

C  trace  computation.   The  Battle  Trace  only  has  meaning  for 

C  battle  periods  which  contain  weapon  effects  related  interactions 

C  between  combatants.) 

BTMAXt ISIDE, IUNIT)  =  0 

ENDIF 
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2  00      CONTINUE 

3  00    CONTINUE 

C Compute  the  Battle  Trace  for  All  Task  Forces  where  threshold 

C  conditions  are  met 

DO  500  ISIDE  =  1,  NUMSIDES 

DO  400  I TASK  =  1,  NUMTASKS 

IF  (  KNUMSHOTS ( ISIDE, ITASK)  .GE.  THRESHOLD  )  THEN 

CALL  BTCOMPUTE( ISIDE, ITASK, PKUNITMAX, PKTASKMAX) 

C Reinitialize  battle  trace  by  setting  accumulator  for  shots  fired 

C        and  recieved  to  zero  after  battle  trace  is  computed  for  a 
C        particular  task  force. 

KNUMSHOTS (ISIDE, ITASK)  =  0 

C Reinitialize  the  accumulators  of  the  influence  of  shots  in  terms 

C        of  SSKP 


C- 


SIPCTPKS( ISIDE, ITASK)   =  0 
SSHOTPKS( ISIDE, ITASK)   =  0 


END  IF 

400      CONTINUE 

500    CONTINUE 

RETURN 
END 
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SUBROUTINE  BTCOMPUTE 


-D.R.  BARR,  DEPT  OF  MATH,  NPS 
J.C.  HOFFMAN,  TRAC-MTRY 
7  APR  92 


SUBROUTINE  BTCOMPUTE (  ISIDE,  ITASK,  PKUNITMAX,  PKTASKMAX  ) 


C- 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c- 


INPUTS: 


ISIDE- 


The  side  of  the  friendly  task  force. 


ITASK-   The  task  force  of  the  friendly  side  for  which 
the  battle  trace  is  to  be  computed. 

PKUNITMAX  -   An  array  with  dimentions  NUMSIDES, NUMTASKS 
which  contains  the  sum  of  the  SSKP  values  which 
constitute  the  contribution  of  all  enemy  units 
to  the  Battle  Trace  measure  of  effectiveness. 

PKTASKMAX  -   An  array  with  dimensions  NUMSIDES, NUMTASKS 
which  contains  the  sum  of  the  SSKP  values 
which  represent  the  capability  of  a  friendly 
task  force  to  contribute  to  the  battle. 

STMETH  -      An  integer  which  designates  the  method 

to  be  used  in  the  computation  of  Battle  Trace, 
rhis  value  is  a  user  input  obtained  from 
CI TV  Janus  Screen  IV  (Form  1119)  . 


INCLUDE 
INCLUDE 
INCLUDE 


'JGLOBE: GLOBAL. FOR' 

' JGLOBE : GLOBUNITS . FOR  ' 

' JGLOBE : GLBTRACE . FOR ' 


DIMENSION 


PKTASKMAX ( NUMS I DES , NUMTASKS ) , 
PKUNITMAX (NUMSIDES , NUMTASKS ) 


Get  Battle  Trace  intermediate  values 

TYPE  *,   '  TIME   ISIDE   ITASK    DELTA  R 

R1DUM  =  SSHOTPKS( ISIDE, ITASK) 

RDUM  =  SIPCTPKS( ISIDE, ITASK) 

B1DUM  =  PKUNITMAX (ISIDE, ITASK) 

BDUM  =  PKTASKMAX (ISIDE, ITASK) 


DELTA  B 


BVAL  ' 


IF  (  BTMETH  . EQ .  1  )  GOTO  100  ! Compute  LOG-TRACE 

IF  (  BTMETH  . EQ .  2  )  GOTO  200  ICompute  Standard  Battle  Trace 

IF  (  BTMETH  . EQ .  3  )  GOTO  300  ICompute  Symmetric  Battle  Trace 

IF  (  BTMETH  . EQ .  4  )  GOTO  300  ICompute  Incrmented  Standard  Trace 

TYPE  *,   'BATTLE  TRACE  NOT  COMPUTED  ***************< 
RETURN 


100 


CONTINUE 


i  ****** 


COMPUTE  LOG-TRACE  ****** 


BTVAL  =  L0G(R1DUM+1)  -  LOG(RDUM+l)  -  L0G(B1DUM+1)  +  LOG(BDUM+l) 
GOTO  90  0 


200 


CONTINUE 


|  ****** 


COMPUTE  STANDARD  BATTLE  TRACE  ****** 


CHECK  FOR  ZERO  ELEMENTS 
IF  (  (R1DUM  .LE.  0)  .OR 
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*  (RDUM  .LE.  0)  .OR. 

*  (B1DUM  .LE.  0)  .OR. 

*  (BDUM  .LE.  0)  )       THEN 
BTVAL  =  1 

ELSE 

BTVAL  =  (R1DUM  /  RDUM)  *  (BDUM  /  B1DUM) 
EIIDIF 
GOTO  900 


300    CONTINUE  ,******  COMPUTE  INCREMENTED  STANDARD  BATTLE  TRACE 

BTVAL  =  ( (R1DUM+.01) / (RDUM+.01) )  *  ( ( BDUM+ . 0 1 ) / ( B1DUM+ . 0 1 ) ) 

******  COMPUTE  INCREMENTED  SYMMETRIC  BATTLE  TRACE  ****** 

IF  (  (BTMETH  . EQ .  3)   .AND.  (BTVAL  . LT .  1.0)  )  THEN 

BTVAL  =  2  -  (1 /BTVAL) 
END  IF 

900    TYPE  *,  CLOCK, ISIDE, ITASK,R1DUM, RDUM, B1DUM, BDUM, BTVAL 

RETURN 

END 
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C SUBROUTINE--BTVIS 

C 


-J.C.HOFFMAN,  TRAC-MTPY 
1  APR  92 


SUBROUTINE  BTVIS  (  IUNIT,  ISIDE,  KOUNT  ) 
C*****   THIS  SUBPROGRAM  SUPPORTS  BATTLE  TRACE  MOE  RESEARCH 


*  C 


C- 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c- 


This  subroutine  builds  a  data  structure  employed  in  the 
computation  of  the  Battle  Trace  MOE. 


CALLING  ROUTINE  --  NRANGE 
INPUTS:     KOUNT  - 


OUTPUTTS 


FUNCTif  N 


USE: 


Number  of  enemy  units  within  visibility 
range 

NMYID  -       List  of  enemy  unit  ID  numbers  within 
visibility  range 

ISIDE  -       Side  of  friendly  unit 

IUNIT  -       Unit  of  friendly  unit 

JSIDE  -       Side  of  enemy  unit  (from  GLOBSRCH) 

Update  of  logical  an  .  .    [NFL  defined  in 
GLBTRACE.FOR 

For  each  observer,  the  appropriate  entries 
are  modified  in  BTINFL  to  indicate  which  enemy 
task  forces  have  at  least  one  unit  which 
is  a  candidate  target  for  the  observer 

BTINFL  values  are  updated  every  call  to  the 
SEARCH  subroutine  in  master  Janus  scheduling 
routine  found  in  RUNJANUS 

BTINFL  contains  logical  data  used  to  compute 
the  values  contained  in  SSHOTPKS  in  the 
subroutine  RELOAD 


-C 
C 
C 
C 
C 
C 
C 
C 

c 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

c 
c 
c 
c 
c 

r 
C 

c 

c 

-c 


INCLUDE 
INCLUDE 
INCLUDE 
INCLUDE 


' JGLOBE : GLOBAL . FOR ' 
' JGLOBE : GLOBUNITS . FOR ' 
' JGLOBE : GLOBSRCH . FOR ' 
' JGLOBE : GLBTRACE . FOR ' 


KTASKFOR 
KOUNT, NMYID 
BTINFL 


200 


DO  2  00  I  =  1,  NUMTASKS 

BTINFL ( ISIDE, IUNIT, I)  =  0 
CONTINUE 


Set  flag  to  zero 


C** 
C 


Set  logical  flag  to  1  for  each  opposing  task  force  which  can  be 
seen  by  the  observer  (IUNIT) 

DO  3  00  J  =  1,  KOUNT 
JUNIT  =  NMYID(J) 
JTASK  =  KTASKFOR (JUNIT, JSIDE) 
BTINFLdSIDE,  IUNIT,  JTASK)  =  1 
3  00    CONTINUE 

RETURN 
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F I LENAME :   GLOBAL . FOR 


INCLUDE  ' JGLOBE:GLBPARAM.FOR' 

INCLUDE  ' JGLOBErGLOBMAIN.FOR' 
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